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The complexation between a semiflexible polyelectrolyte and an oppositely charged macroion 
leads to a multitude of structures ranging from tight complexes with the chain wrapped around the 
macroion to open multileafed rosette-like complexes. Rosette structures, expected to occur for short- 
ranged attractions between the macroion and the chain, have now also been seen in recent Monte 
Carlo simulations with long-range (unscreened) interactions [Akinchina and Linse Macromolecules 
2002, 35, 5183]. The current study provides scaling theories for both cases and shows that rosette 
structures are indeed quite robust against changes in the ionic strength. However, the transition 
from the wrapped to the rosette configuration has a dramatically different characteristics: The 
short-range case leads to a strongly discontinuous transition into a rosette with large leaves whereas 
the long-range case occurs in a continuous fashion. We provide the complete diagram of states for 
both cases. 

PACS numbers: 




I. INTRODUCTION 

Many biological processes involve the complexation be- 
tween a charged macroion and an oppositely charged 
chain. A prominent example is the complexation of DNA 
with histone proteins which is the basic step of DNA 
compaction into chromatin in the cells of animals and 
plants [0. Complexation between macroions and syn- 
thetic polymers is also encountered in many technological 
applications as a means to modify macroion solution be- 
havior; examples are the complexation of polymers with 
charged colloidal particles || and charged micelles [Q . 

There is now a large body of theoretical studies on this 
set of problems, most of them have appeared within the 
last three years |||, ^ |7, ^, ^, 10, 
(cf. also earlier related studiesj 
vary widely with respect to the methods and the level of 
approximations used and also with respect to the physi- 
cal properties assumed for the chain and the macroion - 
the latter usually being niodelled as a charged sphere. 
Some studies g [l| assume the chain to be so 
highly charged that counterion release is the dominant 
mechanism for the sphere-chain complexation, whereas 
most of the other studies assume weakly charged com- 
ponents attracted via standard electrostatics. In some 
cases the ball-chain complexes interact via a short- 
range attraction pO[ |l^, |l^ corresponding to high 
ionic strengths |^|. The investigated systems also vary 
with respect to the chain flexibility: semiflexible chains 
are considered in Refs. §, |, |, 0, 
whereas the other studies assume flexible polymers. Ref- 
erences [^S], ^ were devoted to the diffusion (reposi- 
tioning) of the complexed sphere along the chain. And, 
finally, complexes between a chain and several spheres 
have been considered The space available here 

does not allow us to give a detailed account of the differ- 
ent approaches as well as of the phenomenons predicted 
like overcharging, conformational symmetries of the ad- 



sorbed chain etc. A critical survey will be provided in a 
topical review that is in preparation [plf . 

Sphere-chain complexes have also been investigated in 
several computer simulations l2^, ^ |2j, |2^, 12^, 



30[ |3^. Wallin and Linse studied the effect of 
chain flexibility line charge density ||2^ and sphere 
radius Q on the geometry of a complex between a 
charged sphere and a polyelectrolyte in a Monte Carlo 
simulation that took counterions explicitly into account; 
they also considered the case when there are many chains 
present Chodanowski and StoU investigated 

the adsorption of a flexible chain on a sphere assuming 
Debye-Hiickel interaction. The case of multisphere ad- 
sorption on flexible ||2^ and semiflexible chains was 
studied by Jonsson and Linse. Messina, Holm and Kre- 
mcr demonstrated that in the case of strong elec- 
trostatic coupling it is even possible that a polyelectrolyte 
chain forms a complex with a sphere that carries a charge 
of the same sign - a phenomenon made possible by cor- 
relation effects due to neutralizing counterions. Finally, 
a recent systematic study by Akinchina and Linse 
focused again on the role of chain flexibility on the struc- 
ture of the sphere-chain complex. 

The overall picture emerging from this multitude of 
theoretical and simulation studies is still not very clear. 
Partly this has to be attributed to the fact that there 
are many free parameters determining the properties of 
the sphere-chain complex, especially the length of the 
chain, its linear charge density and persistence length, 
the macroion radius and charge and the screening length 
of the salt solution. This makes it difficult to develop a 
theory that covers the whole range of possible structures. 
The current study is an attempt to give a scaling theory 
that allows an approximate treatment of the chain-sphere 
complex over the whole range of parameters and that 
especially fills the gaps that were left open by the existing 
theories and simulations. It allows to identify the few 
independent scaling parameters in this system and leads 
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to the construction of two-dimensional phase diagrams 
(one for short and one for large screening lengths) that 
cover the whole range of all the other parameters. 

This paper has been induced by a comparison of 
Refs. ^ and |3|. In the former paper we studied 
the complexation behavior of a semiflexible chain and 
a "sticky" sphere. In the context of electrostatics this 
corresponds to high salt concentrations where the Debye 
screening length is much shorter than the sphere radius. 
In that paper we calculated the zero-temperature config- 
urations modelling the polymer by the worm-like chain. 
We found two typical structures that occur in this sys- 
tem. If the sphere is sticky enough, i.e., if the chain 
adsorption energy per length is large, the chain wraps 
around the sphere - as long as there is enough surface 
available. When the adsorption energy is decreased there 
is a point at which the bent chain unwraps in a strongly 
discontinuous fashion. The new structure that emerges 
has multiple point contacts between the sphere and the 
chain and large low-curvature loops connecting them. We 
called this class of structures the rosette configurations. 

On the other hand, the latter work, Ref. presented 
a Monte Carlo study of the complexation of a semiflex- 
ible charged chain with an oppositely charged ball that 
carries the same absolute charge as the chain (isoelec- 
tric complex). No small ions were present so that the 
charged monomers were attracted to the sphere via an 
unscreened long-range 1 / r- interaction. The authors sim- 
ulated systems with different persistence lengths, linear 
charge densities of the chain and ball radii. Depending 
on the choice of parameters they encountered a multi- 
tude of structures - ranging from collapsed structures 
with a "tennisball seam pattern" or solenoid arrange- 
ment of the wrapped chain to open multi-leafed struc- 
tures very much resembling the ones found in Ref. [ p^ . 
Their rosette structures occur for stiffer chains on smaller 
spheres. 

That work demonstrates that the rosette structure is 
quite robust with respect to the range of interaction. At 
the same time it also hints towards major differences be- 
tween these two cases: The short-ranged case predicts 
clearly a strongly discontinuous unwrapping transition 
into large-leafed rosettes when the adsorption energy is 
decreased or the chain stiffness is increased In the 
simulation ||3^, however, the rosette evolved continuously 
with increasing chain stiffness from a tightly wrapped 
state via slightly more open structures with many small 
loops to large-leafed rosettes. 

How can one reconcile these findings? For that purpose 
we reconsider in Section II the chain-sphere complexation 
for the case of short-ranged attraction |l^, formulating 
now the interaction in the language of strongly screened 
electrostatics instead of some short-ranged stickiness. 
Then, in Section III, we give the scaling description in 
the case of weak screening. By focusing on the unwrap- 
ping transition we contrast the two cases in Section IV. 
In the final section we compare our findings to computer 
simulations and experimental observations on the DNA- 



histone complex. 



II. CHARGED ROSETTES AT HIGH IONIC 
STRENGTH 

Consider a polymer chain of length L, radius r and 
persistence length Ip. The chain carries negative charges 
on its backbone - a distance b apart - which leads to a 
linear charge density —e/b. The macroion is modelled as 
a sphere of radius R that carries a positive charge Z . The 
reduced electrostatic potential is assumed to be smaller 
than unity everywhere so that the electrostatic interac- 
tion can be described by standard Debye-Hiickel theory: 
Two elementary charges at distance r interact with the 
potential e(j)/kT = ±lBe~'^''/r with Ib = e^ /eksT (e: 
dielectric constant of the solvent, /csT: thermal energy), 

the Bjerrum length, and k^^ — (SttZbCs)"^^^, the Debye 
screening length for monovalent salt of concentration c^. 
Ib is of the order of 7A in water at room temperature. 

In this section we assume strong screening such that 
the screening length is much smaller than the radius of 
the ball, nR » 1. We will assume throughout the chain 
to be so thin that r is the smallest length scale in our 
system, especially that always nr 1. The adsorption 
energy per length can then be estimated from the Debye- 
Hiickel electrostatic potential close to the surface of the 
sphere which mimics for kR ^ 1 that of a plane with 
a surface charge density a = Z/ (47ri?^). One has then 
e(t)/kT = lBZK~^e~'^^ /R^ where z is the distance from 
the surface (cf., for instance, Ref. |^). This leads to the 
following adsorption energy per length (in miits of fcsT) 



A 



IbZ 
bnR^ 



(1) 



Now the chain will only wrap on the sphere when A is 
large enough, namely so large that it exceeds lp/2R^, 
the bending energy per length |Q . The (free) energy of 
the wrapped chain-sphere complex is then given by 



wrap IpL 



-\L = 
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(2) 



This leads to the prediction of an unwrapping transition 
at X ~ Ip/R^ when the chain stiffness 



IbZ 
bn 



(3) 



is reached. Since the wrapping path of the chain on the 
sphere can be quite intricate (cf. Ref. |^) the local ra- 
dius of curvature is not always precisely R (but of the 
order R) . Thus we dropped from Eq. |^ on all numerical 
prefactors and will also do so in the rest of the paper 
(for similar reasons). Such an unwrapping transition has 
been first discussed by Marky and Manning using 
some unspecified short-range attraction and by Netz and 
Joanny Q for the electrostatic case (in fact, their Eq. (35) 
coincides with Eq. p|). 
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For larger persistence lengths than the one given in Eq. 
^ the chain has to completely unwrap from the sphere. 
This led the authors of Ref. to the prediction of an 
"all-or-nothing" picture: either the chain is wrapped or 
it is unwrapped - with a single adsorption point or even 
completely desorbed. We showed however in a later study 
that the chain can lower its energy considerably by hav- 
ing multiple point contacts to the sphere which leads to 
the rosette structures ||lOj. Let us call fx the energy per 
contact point (in units of kgT). Its value 



XVRi? 



IbZ 



b{KR) 



3/2 



(4) 



follows from the length V Rk~^ of a chain portion around 
a point contact that is located within the distance k^^ 
from the sphere. The free energy of an A^-leafed rosette 
has then approximately the following form 



^rosette 



(5) 



The first term describes the energy required to bend the 
N chain pieces between the point contacts, each of length 
L/N, into leaves with typical curvature L/N, the second 
term account for the TV -I- 1 point contacts. We do not 
account here for the entropy of the chain configurations 
that can be safely neglected for large point contact en- 
ergies /i 3> 1 (see below) . The optimal number of leaves 
follows from the minimization of the (free) energy to be 



ip 



which leads to leaves of size 



eaf 



L 



Ip 



(6) 



(7) 



As already pointed out in Ref. [|lO| each leaf minimizes 
its bending energy by assuming the shape of a so-called 
Yamakawa-Stockmayer (YS) loop, namely a lemniscate- 
shaped loop with an 81-degree apex angle |Q (cf. also 
Ref. ||l6| for an approximate derivation of the YS loop). 
As a cautionary remark we mention that the bending en- 
ergy of YS loops leads to the numerical prefactor 14.04 
for the first term in Eq. || - clearly demonstrating 
that one should not extract from this theory any num- 
bers but just scaling relations. The leaves show negligible 
shape fluctuations for keaf 'ti Ip, i.e. for /i 3> 1. Then 
also N* , Eq. |[ is essentially fixed and the rosette has a 
well-defined leaf number - as implicitely assumed here. 

In Fig. 1 we summarize our results in a "phase di- 
agram" of total chain length L vs. Ip/fi (which is in 
the case of rosettes nothing but the leaf-size, cf. Eq. |^). 
For small values of Ip/fJ. we have wrapped structures, 
for large values rosette conformations. The unwrapping 
transition occurs at Ip/^i ~ R^/kR (cf. Eqs. ^jand^ in a 
strongly discontinuous fashion indicated by a thick ver- 
tical line. The other lines indicate transitions between 



desorption 




FIG. 1: The sphere chain complex at high ionic strength, 
kR 3> 1. Depicted is the "phase" diagram as a function of 
the total length L of the chain and of its persistence length 
Ip divided by the point contact energy fi. The thick vertical 
bar denotes the discontinuous unwrapping transition. To the 
left are the wrapped complexes, to the right the rosettes. 



different ground states, namely rosettes with different 
numbers of leaves. 

Additional features arise from the fact that the sphere 
surface is finite. As we will show now this can lead to the 
formation of tails for the case of wrapped structures and 
to the melting of rosettes. Consider first a wrapped chain 
configuration, i.e., assume some fixed value for Ip/ii < 
R\/kR, and increase L, i.e., go along a vertical line in 
the diagram. Fig. 1. It is clear that at a certain point the 
wrapping has to stop. Locally the wrapped chain forms 
a nearly planar lamellar phase with a distance d between 
neighboring chain sections. This distance follows from 
the competition between the chain-sphere attraction and 
the repulsion between neighboring chain segments [|[ |3^ , 



1 _ i?^ 
~ ba~ bZ 



(8) 



This leads to a wrapping length I ~ R^ /d ~ bZ kso 
which is just the isoelectric wrapping length, the length 
at which the wrapped chain portion just compensates the 
sphere charge. Chains that are shorter than this length, 
L < liso, will be completely wrapped, chains that are 
longer, L > Uso, will have their extra length L — 
dangling from the complex in the form of one or two 
tails. The tail formation ai L — li^o has been indicated 
in Fig. 1 by a horizontal dashed line. We also note that 
if d > in Eq. 0, i.e., if Z < R^n/b (large sphere), the 
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lamellar spacing is reduced to the value as pointed 
out by Netz and Joanny [||. 

Now let us consider th e ro settes. Choose some arbi- 
trary value for Ip/n > Ry/nR and increase L. Then each 
time when the chain length is increased by an amount 
Ip/fi an additional leaf is formed. However, packing con- 
straints imply an upper limit for the number of loops of 
the order iVmax ~ hso/V Rn~^\ this is the maximal num- 
ber of contacts, each excluding an area « dy/RKT^, that 
can be closely packed on the surface of the sphere. For 
a rosette with the maximal number of leaves (in Fig. 1 
we chose arbitrarily iVmax = 5) an increase in L leads 
to an increase of the leaf size, keaf reaches the persis- 
tence length when L ~ A^max^p- At that point melting 
of the rosette takes place, i.e., for larger values of L the 
leaf size distribution is heterogeneous, the YS loops show 
shape fluctuations and also the tails start to grow. En- 
tropic effects are then important. The thermodynamics 
of this melting process is quite intricate and can be cal- 
culated by mapping this problem on an exactly solvable 
one-dimensional many body problem jl^, ^ . In the cur- 
rent study we will not discuss the rosette melting further 
and merely indicate it in Fig. 1 by an arrow pointing 
towards the direction of the L-axis. Also shown in this 
figure is the direction where the chain desorption is to 
be expected. It should occur when the complexation free 
energy of the rosette, Eq. ||, becomes smaller than fc^T, 
i.e. for Ip/fj, > \/lpL. 



III. THE ROSETTE STATE AT LOW IONIC 
STRENGTH 

The rosette configurations discussed in the last section 
occur for chains that are so stiff that wrapping would be 
too costly. It allows a small fraction of the chain to be 
close to the sphere in the form of point contacts. The ma- 
jority of the monomers resides in the loops that do not 
"feel" the presence of the sphere but are needed to con- 
nect the point contacts via low curvature sections. At 
first sight one might thus expect the rosette configura- 
tions to be a special feature for chain-sphere complexes 
with a short ranged mutual attraction. 

That this is not true has already been pointed out in 
the introduction. As observed in the simulations by Ak- 
inchina and Linse |^ rosettes are in fact quite robust 
and occur also in systems with a much larger range of 
interaction. We shall give now a scaling theory that pro- 
vides us with a diagram of states for the chain-sphere 
complex for the case of weak screening, kR <C 1. We will 
first treat the case of short chains L = bN < bZ ~ liso 
where the chain charge is smaller than (or equals) the 
sphere charge. The free energy of the rosette with N 
leaves is then approximately given by 



F, 



rosette 
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knT 



-N 



(9) 



ready discussed after Eq. |^. The second term accounts for 
the attraction between the ball charge Z and the chain 
charge L/b over the typical distance L/N. All other 
contributions to the electrostatics like the monomer- 
monomer repulsion are smaller and neglected in Eq. |^. 
Remarkably this free energy has the same form as the 
free energy |^ of the rosette for short-ranged interaction. 
We just have to identify the point contact energy ^ (for 
the strong screening case) with the single leaf-sphere in- 
teraction 



IrZ 



(10) 



for the unscreened case. The optimal leaf number is thus 
again given by Eq. S a nd the leaf size by Eq. ^ but now 
with fl given by Eq. [lO[ Note that these scaling results re- 
main even true if the leaves grow so large that their outer 
sections do not interact with the sphere due to screening, 
f^heaf > 1- This is so because the electrostatic rosette- 
sphere attraction still scales like —iBZN/b: a fraction 
/ {L/N) (i.e. K~^/b charges) of each of the N leaves 
interacts with the sphere at a typical distance k~^. 

The rosette state competes with the wrapped state. 
We expect that the rosette state is continuously trans- 
formed into the wrapped state when li^af — R] then the 
leaves become so small that they touch with their contour 
the surface of the sphere. Indeed, by setting N = L/R 
in Eq. ^ one finds 



^ wrap 
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(11) 



which can be considered as the free energy of the wrapped 
state: The first term is the bending energy required 
to wrap the chain around the sphere and the second 
accounts for the electrostatic attraction between the 
wrapped chain and the sphere. All other electrostatic 
contributions (as discussed in detail by Nguyen and 
Shklovskii |11 ) are less important and do not occur on 
this level of approximation. On the right hand side of 
Eq. |l^ we arranged the terms in such way that one 
can deduce directly an unwrapping transition, namely 
at Ip/R^ = IsZ/bR which can be rewritten as 



flR: 



IpZR 



(12) 



The first term is the bending energy of N leaves as al- 



Comparing Eqs. || and 111] one might expect that the 
chain unwraps in a strongly discontinuous fashion as dis- 
cussed in the previous section. However, this "unwrap- 
ping point" corresponds just to the point lieaf — R when 
small loops start to grow on the sphere; we thus expect 
the transition to be continuous as pointed out before 
Eq. |ll|. That the unwrapping transition occurs rather 
smooth at low ionic strength and in a strongly discon- 
tinuous way at high ionic strength has been predicted by 
Netz and Joanny however, in that study the authors 
did not allow for rosette structures. 

To complete the picture we finally have to consider 
chains that are longer than the isoelectric length, L > 



liso- We need then at least three terms to capture the 
essential physics of the rosette state: 



N =4 

^ » ma X — ' 
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(13) 



Here we " allow" the monomers to distribute between the 
rosette of length I and a tail of length L — I. The first 
two terms in Eq. |l3| are then as above, Eq. H, the last 
term describes the repulsion of the monomers that con- 
stitute the rosette. The contributions from the tail can 
be neglected. Minimization with respect to I leads to the 
optimal rosette length 



r 



IpN 



h 



(14) 



The free energy |1^ with the optimal length I*, Eq. 14 
minimized for the following number of leaves: 



IS 



N* = 



(15) 



Hence 



this 



level of approximation - I* ~ 
i.e., the rosette monomers just com- 



h^lpN*/lB ^ kso, 

pensate the central ball charge; the rest of the monomers 
extends away from the rosette in one or two tails of total 
length L — kso- Each leaf is of size 

heaf - — (16) 

The rosette disappears for kso/N* = R, i.e., when Eq. [ij 
is fulfilled. It is then replaced by a wrapped structure of 
length liso a-iid one or two tails having the total length 

L liso • 

In Fig. 2 we depict the complete diagram of the sphere- 
stiff chain complexes to be expected in the unscreened 
case. We again plot L vs. Ip/ IJ., the leaf size. When one 
starts in this diagram at a large value of Ip / fi and goes to- 
wards smaller values of (with some arbitrarily fixed value 
L < kso) then all leaves shrink and more and more leaves 
can form. At Ip/fi = R the maximal number of leaves 
(for that given value of L) is reached and at the same time 
the leaves disappear simultaneously in a continuous fash- 
ion. For Ip/fi < R the chain wraps around the sphere. 
For long chains, L > kso, the excess charges are accom- 
modated in tails and all rosettes have the same length 
liso- The borderlines between different rosette ground 
states are then independent of the total length of the 
chain and thus appear as vertical lines. 

Desorption occurs when the free energies, Eqs. ^ and 
p^ , equal the thermal energy fc^T. This point is reached 
when Ip/fi ~ VTpL for short chain, L < kso, and when 
Ip/ fi ~ \/lpkso for long chains, L > kso- We indicate in 
Fig. 2 the direction where desorption occurs by an arrow. 



IV. UNWRAPPING AT HIGH AND LOW 
IONIC STRENGTH 

We take now a closer look at the unwrapping transi- 
tion, contrasting the short- and the long-range case. The 



desorption 




FIG. 2: Diagram of states for the case of low ionic strength, 
K,R <^ 1. The axes are chosen similar to Fig. 1 with fi 
now being the leaf-sphere attraction in the rosette structures. 
The unwrapping is here rather smooth (dashed vertical line). 
Complexes with long chains, L > Uso show tails. 



former case, kR ^ 1, is depicted in Fig. 3(a). As dis- 
cussed in Section II we expect the unwrapping transition 
to occur at A ~ / R^ which leads to Eq. ^. At this point 
the structure jumps in a strongly discontinuous fashion 
into a /an7e-leafed rosette with leaves of size Ip/fi, cf. 
Eq. 1^. Combining Eqs. ^, ^ and ^ we find indeed 



leaf 



(17) 



As discussed in Section II each leaf has the shape of a YS- 
loop iQ with an 81-degrees apex angle. This means that 
neighboring leaves have a relative orientation 9 = 180° — 
81° = 99°, cf. Fig. 3(a). In addition, the leaves have to 
be slightly twisted (like propeller blades) to account for 
the mutual excluded volume. 

The unwrapping at low ionic strength is depicted in 
Fig. 3(b) and goes as follows (we discuss here the case 
L < liso', in the opposite case one has just to replace L 
by liso)- When the chain becomes so stiff that Ip/R^ > 
IbZ/ (bR) the wrapped state is not stable anymore (cf. 
Eq. |ll|). At that point many small leaves {N* ^ L/R 
ones) form simultaneously in a continuous fashion. Their 
size scales as heaf — R, the precise prefactor being not 
accessible to our scaling argument. The number of wind- 
ings around the ball scales as L/R (but note that the 
geometry of the wrapped path can be actually quite intri- 
cate as demonstrated in Ref. Q ) and the typical opening 
angle 7 of each loop at the point of its formation scales as 
(L/R) /N* « 1, again with an unknown numerical value. 
A multi-leafed configuration slightly above the unwrap- 
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place at Ip ~ VlB? /b. The point contact energy is of the 
order {Vl/h) y/RXcc where the so-called Guoy-Chapman 
length Xgc — 1/ (ctIb) is the thickness of the layer of 
condensed counterions around the sphere (which is al- 
ways much smaller than R when there is strong coun- 
terion condensation). The leaf size at the unwrapping 
point is then given by 

Ziea/^ WT^i?»i? (19) 

V Xgc 

which indeed indicates a strongly discontinuous unwrap- 
ping transition. 



FIG. 3: The unwrapping transition (a) at high ionic strength, 
kR ^ 1 and (b) at low ionic strength, kR <^ 1. In the for- 
mer case the transition is discontinuous and large Yamakawa- 
Stockmayer loops are formed, in the ladder case the rosette 
structures evolves continuously. 



ping point is depicted schematically on the right hand 
side of Fig. 3(b). 

Additional insight can be gained by generalizing the 
attractive force between a given chain charge and the 
sphere by a power law —AZ/r" with an arbitrary expo- 
nent a > 0. An integer value a — D~2 with £> = 3, 4, ... 
can be interpreted as a charged chain that adsorbs on an 
oppositely charged D-dimensional ball in D dimensions. 
The electrostatic term for the rosette in Eq. ^ takes then 
the form —AZN"/ (bL"'^) and the one for the wrapped 
state scales as —AZL/bR". Unwrapping takes place at 
Ip ~ AZ/ {bR°'~^) . At this critical value the free energy 
of the rosette Frosette {Ip ^ Ip) has (as a function of N) 
a minimum at N* ~ L/R for a < 2 (Z? < 4), suggest- 
ing a rather smooth unwrapping transition similar to the 
one depicted in Fig. 3(a). This minimum turns into a 
maximum at a = 2 (D = 4). For larger values of a we 
find N* — 0, i.e., the unwrapping transition is strongly 
discontinuous, similar to the short-ranged case discussed 
in the previous section. 

Finally we note that in the case of very highly charged 
chains and spheres the Debye-Hiickel approximation 
breaks down and nonlinear phenomenons like counte- 
rion condensation become important |3^ . The dom- 
inant contribution to the complexation energy is then 
the release of counterions that were condensed on the 
sphere and on the chain before complexation. This is a 
short-ranged interaction and consequently the unwrap- 
ping transition is expected to be discontinuous - even at 
low ionic strength. Following Rcf. 1 14 the wrapped state 
can be written as 



^ wrap 

knT 
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where 51 is the free energy gain (in units of ksT) per ad- 
sorbed chain charge which follows from the entropy gain 
due to counterion release (f2 is a number of order one; 
for details cf. iQ). Unwrapping into the rosette takes 



V. DISCUSSION AND CONCLUSION 

We compare now our results with Monte Carlo simu- 
lations by Akinchina and Linse and then with the 
behavior of a biological chain-macroion complex, the nu- 
cleosome The simulated systems were always at 
the isoelectric point, i.e., the chain length was given by 
L = bZ = ligo- Furthermore, the charges in the simula- 
tion interacted via a non-screened Colombic 1/r poten- 
tial. The simulation results have thus to be compared 
with Fig. 2. Four systems have been considered, each 
having having a fixed set of parameters 6, Z and R but 
with 7 different values of ^p. This means that for each 
case the systems were located on the dashed horizontal 
line at L = liso in Fig. 2. 

In one system (called system II in Ref. [Q) the con- 
tinuous development from a wrapped to the rosette con- 
figurations has been seen most clearly. Example config- 
urations are shown in Fig. 1, system II in that paper 
(that these are representative can be seen by inspecting 
the adsorption probability of monomers as a function of 
the monomer index, cf. Fig. 3 in Q). For /p = 7A the 
chain is wrapped, at Ip = 60A there is already a slight 
indication of very small loops {N = 4 or 5, cf. the small 
oscillations in Fig. 3, systems II, open squares). The 
next system depicted has already a much stiffer chain, 
Ip = 250A, and shows very clearly three leaves, then two 
leaves at Ip = 500A and one leaf for the stiffest chain, 
Ip = lOOOA. In Fig. 2 we have chosen the parameters 
such that iVinax — liso I R equals 4 so that this corre- 
sponds roughly to system II in |^^. To compare with 
the simulations we have to follow the L — ^iso-line in 
Fig. 2: One starts with wrapped structures for Ip/ ^ < R 
and finds then the continuous evolution of rosettes when 
the line Ip/fJ, = R is crossed. The leaves grow at the 
expense of their number (first 4, then 3, 2 leaves), just 
as it has been observed in the simulations. 

In another system (system I) all parameters were kept 
the same except the sphere radius that was now twice as 
large. This system shows 3 small loops when Ip — 250 A 
is reached and one loops for Ip = 500A and lOOOA. In our 
" phase diagram" the transformation R 2R means that 
the unwrapping transition line moves to the right and 
consequently the maximal number of loops goes down by 
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a factor of two: A^max — > which is in satisfac- 

tory agreement with the finding in js^. Going from our 
reference system (II) to another case (system IV) means 
to go to a higher hne charge density of the chain, i.e. 
the transformation b — *■ 5/2 (and thus li^o — * hso/'^)- In 
that case the horizontal hne L — hso has to be moved 
to smaller value which leads again to iVmax A,„ax/2. 
The complex in Ref. js^ shows an unwrapping only for 
Ip > 500A and has already a well pronounced loop at 
Ip = lOOOA. Apparently, the wrapped state is more sta- 
ble in this system than in system I - and the scaling 
argument predicts indeed the critical persistence length 
to scale as 1/b, cf. Eq. The last system (III) finally 
combines a large sphere, R 2R, with a short chain, 
& ^ 6/2. In that case no indication of loops has been ob- 
served. Even at the largest persistence length the chain is 
touching the sphere over some finite section (similarly to 
the geometry discussed by Netz and Joanny g]). In fact, 
we predict A'max Afmax/4 by going from the reference 
system (II) to this case. Altogether, our scaling approach 
is in quite satisfactory agreement with the observations 
made in the Monte Carlo simulations |3^. Interesting 
would be to test our predictions also for chains that are 
shorter or longer than the isoelectric length. Correspond- 
ing simulations are already on the way |^ 

Now we compare our scaling results with an exper- 
imental system, the nucleosome The wrapping of 
DNA around protein spools, octamers of so-called his- 
tones, is the basic step in the condensation of DNA into 
the dense chromatin complex inside the nuclei of plant 
and animal cells. Each repeating unit of this complex, 
wrapped DNA, histone octamer and a stretch of linker 
DNA (connecting to the next such protein spool) is called 
nucleosome. When the linker is digested away the re- 
maining complex, the so-called nucleosome core parti- 
cle, consists of 147 basepairs (« 500A) DNA wrapped in 
l-and-3/4 left-handed superhelical turns around the oc- 
tamer. The core particle which has a radius of 50A 
and a height of '--^ 60A is documented in great detail on 
the basis of high-resolution x-ray analyses . The his- 
tone octamer contains 220 basic side chains (arginine and 
lysine) Q. From these are about 103 located in flexi- 
ble histone tails that dangle of the core par ticle and/or 
are complexed with the wrapped DNA |42|. The rest, 
117 residues, are in the globular part of the octamer, of 
which 31 are exposed to the solvent, the rest being in- 
volved in intra- and interprotein ionic interactions. On 
the other hand, one has 147 bps of DNA wrapped around 
the octamer, each contributing two phosphate groups. 
Hence there are 294 negative charges from the DNA ver- 
sus 220 positive charges of the octamer, i.e., the nucleo- 
somal complex is overcharged by the DNA. 

Yager, McMurray and van Holde |^ characterized the 
stability of the nucleosome core particle as a function of 
the salt concentration (NaCl). Using a variety of ex- 
perimental methods (e.g. gel electrophoresis) they ar- 
rived at the following main conclusions: The core par- 
ticle is stable for ionic strengths ranging from 2mM to 



750mM (this includes physiological relevant salt concen- 
trations around lOOmM). For slighly higher salt con- 
centrations the DNA is partially dissociated; an equilib- 
rium between histones, free DNA and core particles is 
observed. At salt concentrations beyond 1.5A/ the core 
particle is completely dissociated into histone oligomers 
and free DNA. On the other end, for very low salt con- 
centration ^ ImM one finds an "expanded" state of the 
nucleosome. Khrapunov et al. | pl| came via fluorescence 
measurements to similar conclusions: For ionic strengths 
between 5 and GOOmM the core particle is intact. At 
larger ionic strength (« 1.2M) the terminal regions of the 
DNA unwrap and two histone dimers are dissociated and 
at a even larger value (~ 1.5M) the remaining tetramer 
leaves the DNA. Finally, at low salt concentrations one 
encounters an open state: the dimers break their contact 
with the tetramer and the DNA termini unwrap. 

The key features of the behavior of core particle DNA 
(neglecting the substructure of the octamer) were recov- 
ered in a numerical study by Kunze and Netz |^. They 
considered the complexation of a charged, semiflexible 
chain with an oppositely charged sphere, interacting via a 
screened Debye-Hiickel potential. The optimal DNA con- 
figuration was derived numerically from the minimiza- 
tion of the energy. For a reasonable set of parameters 
(comparable to the one of the core particle) they found 
for vanishing ionic strength (k~^ oo) an open, planar 
configuration where only a small fraction of the chain is 
wrapped whereas the two tails (of equal length) are ex- 
tended into roughly opposite directions. This is reminis- 
cent of the open structures reported in the experimental 
studies ^ . Upon addition of salt the structure goes 
from a two- to a one-tail configuration and the chain be- 
gins to wrap more and more. Already well below phys- 
iological ionic strengths the chain is almost completely 
wrapped. The chain stays in this wrapped state up to 
very high salt concentrations. Only then the chain un- 
wraps in a discontinuous fashion when the chain-sphere 
attraction is sufficiently screened. Also these features of 
the complex (a wrapped compact state in a wide range 
around physiological conditions and unwrapping at high 
salt content) reflect the experimental flnding in Refs. 
and ||. 

How are these findings related to our scaling study? 
Let us start at physiological salt concentrations around 
lOOmM. Then ~ 10 A is much smaller than the par- 
ticle size « lOOA and we have the short-range case. Since 
r ~ lOA, i.e., Kr w 1 for physiological conditions, and 
since the binding sites between DNA and the histones are 
quite specific ]40|] , the estimation of A in Eq. ^is not reli- 
able. A can be derived instead experimentally from com- 
petitive protein binding to nucleosomal DNA |Q to be of 
the order (l/5)fcsT/A which is roughly 10 times smaller 
than what Eq. H would predict (but note also that contact 
is made mainly between DNA minor grooves which are 10 
basepairs apart!). This short-range attraction is largely 
balanced by the strong bending contribution {Ip = 500A 
for DNA and hence Ip/ {2R^) ~ 1/6 - l/7fcBT/A). In 
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any case, Eq. ^ predicts an unwrapping for sufficiently 
small values of but the numbers are not reliable. The 
more interesting case is that of the core particle at low 
ionic strength around ImM {k~^ ~ lOOA). At that salt 
concentration 2kR « 1 and we enter the long-range case. 
Apparently, the completely wrapped configuration is not 
stable at this point anymore. This is not surprising since 
the nucleosomal DNA overcharges the protein octamer 
by at least 74 negative charges (if not by ~ 160 charges 
since 86 charged residues are buried inside the octamer). 
We thus expect a considerable part of the terminal DNA 
to unwrap and to be part of one or two tails. In Fig. 2 
this corresponds to the wrapped chain structures with 
tail that are found for small values of Ip/fi, Ip/^^ < R, 
and large values of L, L > kso- 

It would be interesting to redo the experiments with 
DNA pieces that are longer than 147 basepairs. For suffi- 
ciently large salt concentrations there might be then the 
occurrence of rosettes (if there is no interference with 
the partial disintegration of the octamer) . For the other 
limit, it might be appropriate to use the argument for 
highly charged chains and spheres as given at the end of 
the previous section. The linear charge density of DNA 
is very high (2 phosphate groups per base pair, i.e., per 
3.4A). Manning theory indeed predicts that coun- 
terion condensation reduces the linear charge density to 
—b/Ib with Ib = 7 A. Also the charge of the histonc oc- 
tamer is so high that counterion condensation is impor- 
tant. For that case we predicted above the unwrapping to 



take place around Ip = VlR^ /h. However, since the DNA 
persistence length is smaller than E? /b « (50 A) /1.7A 
~ I500A we expect the wrapped state to be stable at low 
ionic strength and no rosettes should occur in this limit. 

Concluding, we have presented a scaling theory for 
the complexation between a charged chain and an oppo- 
sitely charged sphere for the two limits of high and low 
ionic strength. For strong screening one encounters with 
increasing chain stiffness a strongly discontinuous tran- 
sition from a wrapped to an open multi-leafed rosette 
structure as already predicted in Ref. ||l^. In the case 
of weak screening one has again wrapped structures for 
sufficiently flexible chains and rosettes for stiffer chains. 
The unwrapping transition, however, is rather smooth 
in this case. That this is a characteristics of a system 
with long-ranged interaction has been further supported 
by extending this theory to power law attraction with 
arbitrary exponents and to the nonlinear limit of highly 
charged systems. We hope that the presented "phase 
diagrams" for the two cases are helpful for further theo- 
retical and experimental studies. 
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